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Abstract 

We propose a simple model of charge and/or magnetic order formation in sys- 
tems containing both localized and itinerant electrons coupled by the on-site, 
spin-dependent interaction that represents Coulomb repulsion and Hund's rule 
(a generalized Falicov-Kimball model). Ground state properties of the model 
are analyzed on the square lattice on a basis of the phase diagrams that have 
been constructed rigorously, but in a restricted configurational space. For 
intermediate values of the coupling constants there are considerable ranges 
of itinerant electron densities where phases with complex charge and mag- 
netic structures of the localized electrons have lower energy than the simplest 
antiferro- and ferromagnetic ones. A strong tendency towards the antifer- 
romagnetic coupling between spins of localized electrons has been observed 
close to half-filling for any density of localized electrons, including situations 
where the magnetic ions are diluted. For small band fillings the ferromagnetic 
coupling between localized spins is predominant. 
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I. INTRODUCTION 



The purpose of this contribution is to introduce a simple model of charge and/or magnetic 
order formation generated by an on-site, spin-dependent interaction of itinerant electrons 
with localized ones, and to provide a preliminary analysis of that model. The model is 
able to capture many essential aspects of magnetism of both the localized and itinerant 
electrons. It is based on the Falicov-Kimball model (FKM), which was originally proposed 
as a model of metal-insulator transitions in mixed- valence compounds 1 . Later on the FKM 
appeared to be suitable for studying some other phenomena as well, as e.g. a tendency 
for charge- density- wave formation in interacting fermion systems 2-4 . Actually, in most of 
applications studied so far the spinless version of the model was explored. As far as we 
know, spin degrees of freedom were taken into account only in a few papers 5-11 , but in all 
of them the interactions between electrons were spin-independent. 

However, many experiments show that a charge superstructure often occurs together 
with a magnetic order 12-14 . In order to describe both within a single model we proposed 
a generalization of the FKM with an anisotropic, spin-dependent local interaction. These 
ideas have also been applied to the Hubbard model 15 , but that work has been performed 
only in the limit of a large spin coupling constant; the authors named this model the "2-band 
Ising-Hubbard model". Hence we call our model the Ising- Falicov-Kimball model. 

The model Hamiltonian is 

H= E E (* - HdKm)d^ a d n ,a + (Ef - [if) E E fm,r,fm,r, 
<m,n> a =1,l m J?=U 

+ UE E d+ >a d mta f+J m , v (1) 

m <T,V=l,l 

+ J E E {dm,odm,afm,-crfrn,-o d^ n(J d Tn ^ a j^ na j m ^ a ), 

where < m,n > means the nearest neighbor lattice sites m and n, a and n are spin — 1/2 
indices, d m ^ (f m ^) and d^ na (f^ a ) are annihilation and creation operators of itinerant 
(localized) electrons, respectively. The on-site interaction between localized and itinerant 
electrons is represented by two coupling constants: U, which is spin-independent Coulomb- 
type and J, which is spin-dependent Ising-type. The later reflects the Hund's rule force. 
The other parameters of the model are: the hopping amplitude t (we will set it equal to one 
for our energy scale), the energy level of localized electrons Ef and the chemical potentials of 
itinerant fi d and localized [if electrons, respectively. Here Ef and [if are given as independent 
parameters, but physical properties of the system depend only on the difference [if — Ef. 
Taking into consideration that fact, the ground canonical phase diagram was chosen to 
be displayed in the variables [id and [if, as if Ef would be equal to 0. However, it is 
equivalent to another representation with [i and [i — Ef as independent parameters (if one 
puts fi d = [if = n). 

Double occupancy of the localized electrons is forbidden, implying the on-site Coulomb 
repulsion [/// between two / — electrons is infinite. Consequently, at a given site the 
/ — electron occupancy is assumed to be nf = ^/,t + n /a < 1 an d the d — electron occupancy 
to be nd = n d> i+n dt ± < 2. So there are 3 states per site allowed for the f — electrons (nf = 0; 
rifrf — 1 and nf^ = 0; nf^ = and nf^ = 1) and 4 states per site allowed for the d— electrons 
(n d = 0; n diT = 1 and n d>l = 0; n d ,i = and n dd = 1; n d = 2). 
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The first three terms of the Hamiltonian (1) form the multicomponent FKM studied 
recently in the 2D case 10 and in the limit of infinite dimensions 8 ' 9 . The last term of (1) 
describes a simplified interaction that couples spins of an itinerant and localized electron 
occupying the same site. The inclusion of this interaction enables one to describe magnetic 
structures (in addition to charged ones) produced by localized electrons and, at the same 
time, the band magnetism of itinerant electrons. 

All single-ion interactions included in (1) preserve states of the localized electrons, i.e. 
the itinerant electrons traveling through the lattice change neither occupation numbers nor 
spins of the localized ones. In other words, [H, f^fi V ] = for all % and 77, so the local 
occupation number is unchanged. This is a characteristic feature of the Falicov-Kimball 
based models, which makes them tractable in a controllable way. 

The localized electrons play the role of an external, charge and spin dependent potential 
for the itinerant electrons. This external potential is "adjusted" by annealing, so the total 
energy of the system attains its minimum. So there is a feedback between the subsystems 
of localized and itinerant electrons, and this is the feedback that is responsible for the long- 
range ordered arrangements of the localized ones, and consequently for the formation of 
various charge and/or spin distributions in low temperatures. 

On the other hand, if the total magnetization of the localized electrons is non-zero, the 
potential experienced by the itinerant spin-up electrons differs from that of spin-down ones. 
As a result, the spins of the itinerant electrons become partially polarized, so we also have 
unsaturated band magnetism. 

So far it was a common practice to investigate the isotropic, Heisenberg-type interac- 
tion between spins of localized and itinerant electrons (e.g. within a framework of the s-d 
model 16 ). Here, instead, we propose an Ising-type coupling between the spins (see also 
Ref. 15 ). An advantage of the latter approach is that the model can be treated rigorously. 
However, one can also provide a plausible physical background that justifies the assumption 
that the spin-flip processes generated by S + S~ and S~S + operators, which are character- 
istic for the Heisenberg-type interaction, can be neglected in a first approximation. The 
reasoning comes from a simple notice that magnetic structures observed in many materials 
are stable. So, presumably, electrons moving through a crystal preserve their spins over 
many lattice sites. Besides, the Ising type coupling between itinerant electrons was already 
studied in the framework of the s — d model 17 . 

Here we show that the model (1) captures driving mechanismes of formation of stripe 
phases and other charge and/or magnetically ordered superstructures. The class of materials 
fitting this picture is quite big and encompasses conductors that are magnetically ordered 
in low temperatures, but with no sign of the Kondo effect (as the latter implies the spin-flip 
process). There is a huge number of such systems, including many compounds of lanthanides 
and actinides, where various complicated magnetic structures were detected 18 . However, the 
model we propose is primarily aimed to desrcibe systems that display both a charge and 
magnetic order. In particular, we expect that dopped systems with a large crystal field effect 
would be the best candidates. The dopped systems, as they are usually composed of ions 
with different occupancies of localized electrons (a mixed-valency regime) and because they 
usually contain band electrons. This category includes materials, where stripe phases were 
recently detected as e.g. Lai,QNdoASr x Cu0^ 19 , Y Bo,2Cu30q +x 12 , Bi 2 Sr2CaCu20s+s 14 or 
La,i. 5 Sr . 5 NiO4 20 . The behaviour of systems like those was already analyzed theoretically 
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on a basis of the extended Hubbard model 21 ' 22 and the t — J model 23 . 

On the other hand, a strong crystal field makes a flip of magnetic moment of an ion 
difficult. Consequently, a process of a simultaneous flip of spins of a localized and an 
itinerant electron becomes rare. This is another justification of taking into account only the 
Ising-type, instead of isotropic Heisenberg-type coupling in the model. 

A major interaction omitted in the Hamiltonian (1) is the Hubbard-type interaction be- 
tween spin-up and spin-down electrons. An inclusion of that term would make it intractable 
rigorously for arbitrary values of the coupling constants, but once again, one can imagine 
a simple justification for the omission. Namely, it is assumed here that the strength of the 
on-site interactions between two particles abides by the following hierarchy: a) the largest 
one - when the both particles are localized, b) an intermediate one - if one of the particles is 
localized and the other is itinerant, c) a negligible one - when the both particles are itiner- 
ant. This hierarchy may be summarized by a rather natural rule: the longer time particles 
occupy the same site, the more important becomes interaction between them. We point out 
here that the cases a) and b) are treated exactly within the model we propose. Obviously, 
there are no obstacles to include various interactions neglected in the Hamiltonian (1) as 
perturbations. On the other hand, the simplicity of the model (1) makes it attractive for 
studies using various techniques, both analytical and numerical. 

An influence of the Hubbard interaction Udd between spin-up and spin-down itinerant 
electrons on an arrangement of the localized electrons is not known in a general case, but 
for Udd small we don't expect any dramatic changes of the phase diagram. However, if Udd is 
large, the situation may be different. For example, for half-filling (one itinerant electron per 
site) the itinerant electrons are ordered aniferromagnetically in the large Udd limit and they 
impose the same simple type of order to the localized electrons. Then, if there is one localized 
electron per site, the configuration of the localized electrons will be antiferromagnetic, i.e. 
the same as in the case studied in the current paper. If, however, the density of localized 
electrons is equal to 1/2, the ground-state configuration D3 (see Fig. 2c) of the model (1) 
will have higher energy than the phase corresponding to the configuration b displayed in 
Fig. 2a. So, if the Hubbard interaction is taken into account transformations of ground- 
state configurations of localized electrons are expected for some pairs of densities of electrons 
{pd,Pf) but not for the whole range of their values. 

In the current paper we study the ground-state phase diagrams of the model (1) on 
the square lattice. However, the method we use can be applied to various types of one- 
two- or three — dimensional lattices. The studies are based on the method of restricted 
phase diagrams 3 constructed in the grand canonical ensemble (in the plane (pd, Pf) - see 
Fig. la) and then translated into a canonical ensemble diagram (in the plane (pd, Pf) - see 
Fig. lb). Working within the framework of a grand canonical ensemble first assures the 
thermodynamic stability conditions are fulfilled (see Ref. 24 ). 

Another important reason for using the chemical potentials pd and pf in the model 
(1) is a possibility of changing and adjusting the electron's occupation numbers. Then, 
although the model does not contain a hybridization term, suitable changes of the chemical 
potentials or the position of the Ef level can produce appropriate changes of the occupation 
numbers. In particular, fixing a total number of electrons enables us to study classical 
intermediate valence states, where some localized electrons change their occupancy to move 
into the delocalized states, conserving the total number of electrons. Of course, a change of 
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the chemical potentials or Ef cannot evoke a change of the electron states but only their 
occupation numbers. 

Previous work on the spinless FKM has revealed a rich structure of the phase diagrams 3 ' 4 . 
A number of charge ordered superstructures has been found, including stripe and nonstripe 
phases. However, so far little is known about magnetic structures accompanying those of 
charge ones, except for the limit of large J, where the ferromagnetic order was proven to 
exist for pd < pf < l 15 . As far as we know, the present study represents the first analysis of 
both charge and magnetic ordering in the framework of a generalized Falicov-Kimball based 
model for finite U and J values. 

In the next section we explain how the calculations were carried out. The restricted 
phase diagrams are presented and described in the section III. The last section contains 
summary and conclusions. 



Here we consider all possible configurations of the localized / — electrons (including 
their spins), for which the number of sites per unit cell N is less or equal to iV c = 4. 
Taking into account all allowable sizes and shapes of the unit cells, as well as all relevant 
translational vectors, we found 47 distinct configurations. Some of them, representing phases 
with the smallest periods (up to 2 sites in an unit cell) and those with higher periods that 
are analyzed in this paper, are displayed in Fig. 2. The presentation of the configurations is 
chosen in such a way that make an easy observation of their transformations, starting from 
the ferromagnetic phase (F) and ending either at the antiferromagnetic one (AF) (Confs. 
1-4), or at the "empty" one (E) (Confs. D1-D5), in accordance with the diagrams given in 



For each periodic configuration in our trial set we performed the Fourier transformation 
of the Hamiltonian (1) and determined the electronic band structure for the conduction 
electrons. In other words, we solved the eigenvalue problem and found the eigenvalues E ver k, 
with branch index v = 1,2,..., Nq, spin index a and the Bloch wavevector k = (k x , k y ) (for 
more details see Refs. 3 ' 4 ). This required us to diagonalize up to 4 x 4 matrices and resulted 
in analytical formulae for at most 4 different energy bands, separately for spin-up and spin- 
down electrons. (It is the main reason why the maximum size of the unit cell is limited to 
N c — 4 only, as the analytical formulae for phases with larger unit cells are not known in the 
general case.) In the simplest cases, related to 6 configurations with unit cells containing 1 
or 2 sites (Confs. E, F, AF, a, b - see Fig. 2a and Conf. 3 - see Fig. 2b), the eigenvalues 
are given in the Appendix. For the remaining configurations, with unit cells containing 3 or 
4 sites, the analytical formulae are quite long and for that reason are not displayed here. 

Having exact formulae for the energy spectra E va \. all quantities of interest can be cal- 
culated from the densities of states of spin-up and spin-down electrons 



II. METHOD OF CALCULATION 



Fig. 1. 




(2) 



In particular, the electron densities 




(3) 
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and the total electronic energy per site 

EtotiHd) = / E{Z } {E) + Z { {E))dE (4) 

J — oo 

as a function of the chemical potential (the Fermi energy) p d . The energy per site in the 
ground canonical ensemble (the Gibbs thermodynamical potential) is defined as 

E gc = E tot - p d (p d] + p di ) - /i/(p /T + Pfi). (5) 

For the simplest configurations E, F and AF the quantities can be obtained analytically, but 
the expressions are rather complicated (they are given in terms of special functions). This 
is way in our calculations we employ a Brillouin zone grid of 100 x 100 momentum points 
for each bandstructure. At each discrete momentum point of the Brillouin zone we get at 
most 4 eigenvalues. Then the eigenvalues of the bandstructure are summed to determine the 
ground-state energy for each number of conduction electrons. The Gibbs thermodynamical 
potential is calculated for all possible values of the chemical potentials of the conduction 
and localized electrons through the formula 

Egc({Wi}) = E( 1Q0 Q 0iV 12 e J*( W) - VdPdo - PfPfa), (6) 

where Wi = or j or j characterizes the localized electron state at site i, and the symbol 
€j a ({wi}) denotes the energy eigenvalues of the bandstructure for the given configuration of 
localized electrons {wi}. The ground canonical phase diagram is constructed in the plane 
of (p d ,pf). We directly compare the ground state energies of all phases from the trial set, 
and select the lowest one. Finally, the grand-canonical diagram is translated to a canonical 
phase diagram for arbitrary densities pf = P/t + P/| and p d = p d \ + p d \ of the /— and 
d — electrons, respectively. This procedure assures thermodynamical stability of all phases 
(both periodic and their mixtures) present in the resulting canonical phase diagrams 24 . 

III. RESULTS AND DISCUSSION 

The grand canonical phase diagrams plot phase boundaries, which divide the plane of 
chemical potentials into domains, in each of which a single localized electron configuration 
forms the ground state. The general structure of the diagrams (Fig. la) is similar to that 
found for the spinless FKM 3 . It is relatively simple for large U and more complex for small U. 
However, since the spin degeneracy is now lifted, new domains related to various magnetic 
arrangements having the same distribution of localized particles appear. In particular, the 
upper region corresponding to the full configuration in the spinless case is now split into 
a set of phases, ranging from F, through ferri- and various complex antiferromagnetic to 
the simplest AF phase; see Figs. 1, 2a,b. The lower region corresponds to the empty 
configuration (E - see Fig. 2a), and the largest of the remaining regions is the diluted 
antiferromagnetic phase (D3 - see Fig. 2c). 

In the case of U = 8, J = 0.5 the central portion of the diagram is divided into diagonal 
stripes in which various periodic configurations of localized electrons are ground states. The 
sequence of the localized electron densities, reading from the left to right at fixed pf is 3/4, 
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2/3, 1/2, 1/3 and 1/4 (Confs. D1-D5, see Fig. 2c), the p f = 1/2 phase (D3) being the 
largest region. 

In the corresponding canonical diagram (in the plane (p d ,pf) ) given in Fig. lb all these 
phases are situated along the line pj + p d /2 = 1. Movement along that line (starting from 
F) corresponds to a dilution of the ferromagnetic state by a gradual replacement of the 
localized electrons by itinerant ones in such a way, that each localized electron is replaced 
by two itinerant ones. This is different from the spinless case, where periodic phases are 
located along the line pf + p d — 1. The difference is related to the fact, that now both spin 
up and spin down itinerant electrons can independently occupy each site, so the total density 
p d can take any value between and 2. 

The analysis of the canonical phase diagram displayed in Fig. lb shows a strong tendency 
to antiferromagnetic order in the region of densities p d around 1 (half-filling) and to the 
ferromagnetic order outside that region. Indeed, both the phases D2 and D3 (see Fig. 2c) 
situated along the line pj + p d /2 = 1 and those represented by the configurations no. 3, 4 
and AF (see Figs. 2a,b) situated along the line pf — 1 are ordered antiferromagnetically. 

For intermediate values of the density p d and pf — 1 both antiferro- and ferrimagnetically 
ordered phases are found (Conf. 1 and 2 - see Fig. 2b). These phases represent consecutive 
stages of a reconstruction process when F transforms into the simplest AF with an increase 
of pd- 

On the other hand, only ferro- (F, Dl, D4, D5 - see Figs. 2a, c) and antiferro- (D2, 
D3 - see Fig. 2c), but not ferrimagnetically ordered phases were found along the line 
Pf + Pd/2 = 1. Moving along that line (in Fig. lb) from the left-upper corner (p d = 0,p/ = 1) 
to right-lower one (p d = 2,pf = 0) one can notice that the process of dissolution of the full 
ferromagnetic phase overlaps with the tendency to antiferromagnetism inside a region around 
half-filling. 

The picture becomes more complicated as U is reduced. An example is the case of U = 2 
and J = 0.5 (not displayed here), where complex structures of the diagrams are found both 
in the plane of chemical potentials (pd,Pf) and densities (pd,Pf). First of all, more phases 
appear in the diagrams (20, versus 12 in the case of U = 8) and the domains occupied by 
the phases in the central region are no longer diagonal stripes parallel to each other, but 
have a less regular structure. 

The phases from the central region are distributed in different sectors of the plane (p d , pf) 
in the corresponding canonical phase diagram. Only three periodic phases are found along 
the line pf + p d /2 = 1: D3, D5 and Dla - see Figs. 2c, d. Two of them (D3 and D5) 
already appeared in the U = 8 case, but the ferrimagnetically polarized Dla phase replaces 
the ferromagnetic Dl (compare Figs. 2c, d). This fact illustrates a possibility of phase 
transitions from one magnetic structure to another with a change of U. 

Another characteristic feature of the periodic phases found in the U = 2 diagram is their 
stability with respect to finite intervals of the density p d . This means that the Fermi level 
of the corresponding phase lies inside its energy band, so those phases are conducting. This 
is just the opposite to the situation found in the U — 8 diagram (Fig. lb), where the phases 
D1-D5 are ground states for fixed values of itinerant electron densities p d , corresponding to 
insulating states (the Fermi levels are situated inside their energy gaps). 

Presumably, if we increase the maximum size of unit cells N c , then phases with larger 
periods will enter the phase diagram in such a way, that intrevals of the density p d where 
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low-period phases are stable will decrease. It is not clear from our studies if the intervals 
shrink to single values or to finite intervals of the density in the full phase diagram. Perhaps 
the first scenario occurs for some phases and the other for the rest. 

A common feature of the canonical diagrams for U = 8 and U — 2 is that in both 
cases the fully occupied phases with pf — 1 (Confs. F, 1-4, AF) are the same and have 
identical positions. It appears that a distribution of these phases is symmetric with respect 
to half-filling pf — 1 and it does not depend on U, but merely on J. 

Phases with p/ = 1 that appear in the diagrams shown in Fig. 1 seem to be stable 
within finite intervals of p d . This means that their ground states can be conducting. In 
particular, antiferromagnetic conductors appear to be possible ground states for certain 
electron concentrations off of half-filling. However, this is probably an artifact due to the 
restriction imposed by the maximum period of configuration in the trial set. If one takes 
into account configurations with larger periods then, presumably, some of them will enter 
the phase diagrams in such a way, that more and more phases will be stable merely within 
energy intervals lying inside one of their energy gaps. Therefore, one expects that the full 
phase diagram will contain many periodic phases in the insulating states, and only a part of 
them, especially those having gapless energy spectra, as e.g. those represented by the Conf. 
E and F, will be conducting. 

A similar situation was observed for the spinless FKM, where large period phases oc- 
cupy regions of the phase diagram located outside of the energy gaps of phases having low 
periods 2,3 . 

IV. SUMMARY AND CONCLUSIONS 

A simple model capable of describing both charge and magnetic structures of localized 
electrons, as well as itinerant band magnetism, was introduced and investigated on a square 
lattice. Restricted ground state phase diagrams were constructed for intermediate values 
of the coupling parameters. Various types of charge and magnetically ordered phases were 
detected for a range of band fillings, illustrating e.g. consecutive stages of transformation of 
F to AF with an increase of the band filling. 

It is remarkable, that the results presented here are consistent with those found for the 
Hubbard model in the large U limit 25 and also with those obtained in Ref. 15 , where - in 
turn - the limit of large J was investigated, as in all these cases one gets the simplest AF 
at half-filling and F far away from that limit. However, according to Ref. 15 and Ref. 25 the 
ferromagnetic ground state extends to all band fillings but half-filling, which is consistent 
with the Nagaoka theorem 26 . This, of course, is a consequence of investigation of the large 
U or J limit. 

Here, instead, we provide results that give an opportunity to determine the upper limit 
for the density of itinerant electrons below which F can be stable for finite U and J. Indeed, 
a direct comparison of energies of phases from the restricted set (they were calculated rigor- 
ously by using the exact analytical formulae for their energy spectra) shows that F can be 
stable for p d < 0.131 (or p d > 2 - 0.131) if J = 0.5 and for p d < 0.233 (or p d > 2 - 0.233) 
if J = 1. These results complete those given in Ref. 15 , where no critical value for p d below 
which F is stable was reported. From the theorems provided in Ref. 15 it may be merely 
concluded, that F is stable for any < p d < pf if pj = 1 and J is large enough. 
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Since no restrictions were here imposed on values of U or J, it was not a great surprise 
that other ordered phases than just the simplest AF and F have been detected for interme- 
diate electron concentrations. In particular, various magnetically ordered phases were found 
for pf = 1 (when charge is distributed uniformly) and intermediate U and J, illustrating a 
transformation of F to AF with an increase of band filling. 

But what are the physical reasons and what are the driving forces that lead to the complex 
charge and/or magnetic arrangements? If one attempts to describe such behavior in terms 
of two-body forces, then interactions between more distantt than the nearest neighboring 
lattice sites needs to be exploited. The double exchange 27 and RKKY interactions, derived 
from the correlated electron models by using perturbative methods, are prominent examples 
of such approaches. 

On a basis of the present studies we conclude that the phenomenon of charge and/or 
magnetic ordering can be explained by a mutual adjustment of the distributions of the 
charges and spins of the itinerant and localized electrons. This idea is not new, since it 
constitutes a background for the famous density functional theory. The difference is that 
we already start from an effective Hamiltonian (represented here by (1)) and consider many 
lattice sites instead of focusing on the details of a charge distribution around a single site, 
as is commonly practiced in the ab initio calculations. This allows us to take into account 
the kinetic energy of the itinerant electrons and to extract the most essential information 
about the system under investigation. 

Indeed, one needs to investigate large enough areas accessible for itinerant electrons, as 
their total energy depends on the distribution of the localized ones over many lattice sites. If 
the latter forms a charge and/or magnetically ordered structure, then the mean values of the 
charge density distribution and/or the spin polarization of the moving electrons also adopts a 
suitable structure. It is clear that in the simplest case of a ferromagnetic metal, the densities 
of itinerant spin-up and -down electrons spread out uniformly over the whole lattice, but 
their values on each site are different, i.e. the moving electrons are polarized. But in the 
general case the distributions of itinerant spin-up and -down electrons are inhomogeneous 
and they may differ one from another. 

Such an inhomogeneous distribution means that the itinerant electrons have some more 
preferable routes of travelling though the crystal (the routes may be different for the spin-up 
and -down electrons). So one can notice here something like a traffic self-regulation, where 
the system of localized electrons orders in such a way that the itinerant electrons would have 
as much freedom to move as possible. This of course is governed by the quantum mechanical 
laws through a minimization procedure of the total energy, which depends not only on the 
coupling constants, but also on the densities of the electrons. 

It is interesting that many of the stable structures are axial (the lines of equivalent sites 
are parallel to the lattice axis) or diagonal (the corresponding lines are oriented along the 
(1,1) direction) stripes. Axial stripes are predominant for rather small band fillings, i.e. 
for small densities of itinerant electrons, and the diagonal stripes are found close to the 
half-filling. The same picture was already observed for the spinless FKM 4 . 

This stripe-type ordering means that the itinerant electrons prefer to move along simple, 
one-dimensional channels. Has this observation anything to do with, or can it shed light on 
a mechanism of the high — T c superconductivity? It is too early to answer this question, but 
the fact is that the stripe ordered phases have been detected in many high — T c materials. 
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So further studies along the lines indicated here will be interesting. 

Another point is that the detection of charge and magnetic phases other than F or AF 
provides an opportunity to describe various structures observed in many systems. Indeed, 
for many years various complex magnetic structures were observed in materials containing 
transition metals, rare-earths or actinides. Recently, thanks to new precise experimental 
techniques, stable charge superstructures have been also found in a number of compounds 
(see e.g. 12 and the references given there). The model proposed here is capable to describe 
such structures in a simplified way but, of course, this work is merely the first step towards a 
complete analysis of the apparently very complicated processes occurring in real materials. 

It is worthwhile to notice, that even though only the restricted phase diagrams were 
analyzed here, the conclusions are expected to also hold for the complete phase diagrams. 
This conjecture comes from the fact that the reported results are consistent with those 
obtained exactly in limiting cases and from a comparison with the data found previously for 
the spinless FKM, where an increase in the size of the allowed unit cells does not produce 
significant qualitative change in the phase diagram. 

In closing we admit that it is clear that the simple effective model presented in this paper 
cannot describe many interesting phenomena observed in solids, like the Kondo effect or the 
superconductivity. However, since it is already a non-trivial model that can describe some 
of the phenomena, it can serve as a reference system for studies of more elaborate models 
of correlated electron systems. 
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Below are given exact formulae for the energy spectra Ej^ of the spin-up itinerant 
electrons in the simplest cases of period 1 and 2 configurations for the localized electrons 
(Conf. E, F, AF, 3, a, b - see Figs. 2a,b). The corresponding spectra of the spin-down 
electrons can be obtained by taking the opposite signs of parameter J. 



APPENDIX A 



Ei k f(k x , k y ) = 2(cosk x + cosky) 
E[ k JU, J, k x , k y ) = U — J + 2(cosk x + cosk y ) 



(Al) 
(A2) 

(A3) 
(A4) 



Ef£(U, J, k x , k y ) = UT + 4:(cosk x + cosky) 2 



Ei k ^(U, J, k x , k y ) = U + 2cosk y =F \J J 2 + A(cosk x ) 2 




(A5) 




(A6) 
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FIGURES 



FIG. 1. The restricted ground-state phase diagrams for U = 8 and J = 0.5: (a) the ground 
canonical phase diagram, i.e. in the plane of chemical potentials (/x^, Pf), (b) the canonical 
phase diagram, i.e. in the plane of densities (pd, Pf). F, AF and E refer to the ferromagnetic, 
antiferromagnetic and empty (i.e. without localized f-electrons) phases, respectively. Periodic 
phases with pj = 1 are denoted by the numbers 1-4 (see Fig. 2b), and those with pf + pa/2 = 1 
by the symbols D1-D5 (see Fig. 2c). (The domains D2 and D4 are very thin, therefore their 
boundaries are so close to each other that they have the appearence of a single thick line.) In the 
case (a), the phases occupy finite size domains and mixtures refer to the lines that separate the 
domains. In the case (b), the phases are ground states only on the bold straight-line segments or at 
single points. Outside these straight line segments or points there are mixtures of periodic phases 
that have lower energy than any periodic phase taken from the restricted set. The small vertical 
straight line segments crossing the pf = 1 line mark limits of F. The diagonal line pf + = 1 
is only a visual guide. 

FIG. 2. Pictures of configurations of the localized electrons representing phases discussed in the 
paper. The f (j) refer to sites occupied by a localized electron with spin up (down) respectively, 
and the dots denote sites not occupied by the localized electrons. The shaded region in the lower 
left corner shows the unit cell, and line segments show the translation vectors that are used to tile 
the two-dimensional plane. Panel (a) shows the simplest configurations of the localized electrons 
containing 1 or 2 sites per unit cell. Three of them (E - empty, F - ferromagnetic and AF - 
antiferromagnetic) appear in the phase diagrams and the last two (a and b) are not present in the 
diagrams. Panel (b), the configurations 1-4 with pf = 1, which are found on the restricted phase 
diagrams for U = 8 and J = 0.5 between F and AF. Panel (c), the configurations D1-D5 with 
Pf < 1, which are found in the restricted phase diagrams for U = 8 and J = 0.5 between the F and 
E phases; on the canonical phase diagram (Fig. lb) they appear along the diagonal pf + pd/2 = 1- 
Panel (d), the configurations with pf = 1 (Confs. la, 5), that are found between the F and AF 
phases in the restricted phase diagrams for U = 8 and J = 1.0, but are not found on the diagrams 
for U = 8 and J = 0.5: Conf. la replaces Conf. 1, whereas Conf. 5 appears between Conf. 4 and 
AF. The configuration Dla with pf = 0.75 is found on the restricted phase diagrams for U = 2 
and J = 0.5 instead of Dl (see Fig. 2b). 
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